Double Quantum Dots in Carbon Nanotubes 



J. von Stecher,* B. Wunsch,* M. Lukin,* E. Dernier,* and A. M. Rey* 
* JILA, University of Colorado and National Institute of Standards and Technology, Boulder, Colorado 80309-0440, and 
' Physics Department, Harvard University, Cambridge, Massachusetts, 20138, 
(Dated: August 31, 2010) 

We study the two-electron eigenspectrum of a carbon-nanotube double quantum dot with spin- 
orbit coupling. Exact calculation are combined with a simple model to provide an intuitive and 
accurate description of single-particle and interaction effects. For symmetric dots and weak magnetic 
fields, the two-electron ground state is antisymmetric in the spin- valley degree of freedom and is 
not a pure spin-singlet state. When double occupation of one dot is favored by increasing the 
detuning between the dots, the Coulomb interaction causes strong correlation effects realized by 
higher orbital-level mixing. Changes in the double-dot configuration affect the relative strength of 
the electron-electron interactions and can lead to different ground state transitions. In particular, 
they can favor a ferromagnetic ground state both in spin and valley degrees of freedom. The 
strong suppression of the energy gap can cause the disappearance of the Pauli blockade in transport 
experiments and thereby can also limit the stability of spin-qubits in quantum information proposals. 
Our analysis is generalized to an array of coupled dots which is expected to exhibit rich many-body 
behavior. 

PACS numbers: 73.63.Fg,73.23.-b,73.22.-f 



I. INTRODUCTION 

Experiments on few-electron double quantum dots al- 
low the measurement and manipulation of the spin de- 
gree of freedom of the confined electrons-M Such con- 
trol is at the heart of semiconductor-based spintronics^' 
and quantum-information proposals^'. Recently, sub- 
stantial experimental efforts have been focussed on con- 
trolling electrons in carbon nanotube double quantum 
dots, and many o f th e capabilities previously achieved in 
GaAs double dots^HH are starting to be reproducecPEJ] 
These include the ability to start from an empty dou- 
ble dot and systematically fill it with electrons. Since 
12 C has no nuclear spin, carbon based nanostructures 
are expected to reduce hyperfinc induced decoherence as 
compared with GaAa^. 

Furthermore, carbon based materials exhibit richer 
physics than GaAs semiconducting materials because of 
the additional valley degree of freedonPHUl. 

In princi- 
ple, the spin and valley degrees of freedom could lead to 
a SU(4) symmetry at zero magnetic field instead of the 
standard SU(2) symmetry in conventional semiconduc- 
tors (see e.g. Ref . [TBI and references therein). However, 
it has been recently demonstrated^ that the enhance- 
ment of the spin-orbit splitting in small-radius nanotubes 
breaks the four-fold degeneracy of the single-electron 
ground state into a two-fold degeneracy. 

In this work, we study how spin-orbit coupling and 
electron-electron interaction effects are manifested in the 
two-electron spectrum and transport properties of a car- 
bon nanotube double dot. This represents an extension 
of previous studies on few-electron physics in a single 
carbon nanotube dotP^Sl. We find that, despite of spin- 
orbit coupling and the existence of an additional valley 
degree of freedom, the two-electron eigenstates can be 
separated in an orbital part and a spin-valley part that 



are, to a very good approximation, independent of each 
other. The spin-valley part can be grouped in six an- 
tisymmetric and ten symmetric spin-valley eigenstates 
which we refer to as multiplets. 

The separation of spin-valley degrees of freedom sig- 
nificantly simplifies the description of the systems and 
allows us to draw analogies with standard GaAs double 
dots. Our main results can be summarized as follows: 
(a) For dots at zero magnetic field and no detuning, each 
dot is populated by a single electron and tunneling is sup- 
pressed because of Coulomb interactions. Thus, interdot 
coupling only occurs virtually via superexchange interac- 
tions that determine the ground state symmetry. In this 
regime, we find a spin- valley antisymmetric ground state, 
that does not have a well defined-spin due to the spin- 
orbit coupling, (b) For large detuning between the dots, 
double occupation of the same dot becomes favorable. 
In this regime, Coulomb interactions can mix higher or- 
bitals in the two-electron ground state^"^. This ad- 
mixture significantly reduces the energy spacing between 
the multiplets (which, for weakly interacting electrons, is 
determined by the orbital level spacing). The interplay 
between spin-orbit coupling and interaction then leads to 
a ferromagnetic ground state above a small critical mag- 
netic field, since the Zeeman terms overcome the strongly 
suppressed splitting between effective singlet an triplet, 
c) The reduction of the energy gap between orbitally sym- 
metric and antisymmetric states caused by the Coulomb 
interaction affects transport properties through the dot 
and can lead to the disappearance of the so-called Pauli- 
blockade (suppression of current through the double dot 
due to the Pauli exclusion principle) and might explain 
the absence of Pauli blockade reported in Ref. [TT] The 
absence of spin-blockade might affect the performance of 
quantum information proposal which use spin-qubits in 
double dots, since gate operation in those proposals is 
based on the spin-blockade mechanism. The disappear- 
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ance of Pauli blockade might be prevented by reducing 
Coulomb correlations by either working with short dots, 
or by covering the nanotube by large dielectrics^. 

This paper is organized as follows. In the section II, we 
introduce the microscopic model for the double dot and 
analyze the non-interacting predictions taking into ac- 
count a magnetic field parallel to the nanotube axis, spin- 
orbit couplings and detuning between the dots. Then, we 
construct a simple two-electron model that captures the 
interaction effects. This model is then compared with 
solutions of an exact many-band Hamiltonian using lo- 
calized single-particle orbitals. In Section |Hl} we discuss 
the energy spectrum of two interacting electrons in a dou- 
ble dot in three different detuning regimes corresponding 
to (i) a symmetric double dot, with one electron in each 
dot, (ii) strong detuning, with both electrons in the same 
dot, and (iii) at the crossover between both regimes. We 
then analyze the transport properties of the double dot. 
Finally, we discuss how to extend our low energy behavior 
analysis to serially coupled quantum dots. In section IV, 
we present the conclusion. Technical details on the cal- 
culation of Coulomb matrix elements and the derivation 
of the rate equations used for transport are presented in 
Appendices [A] and [B] 



II. MODEL 

Single-wall carbon nanotubes are formed by a single 
layer of graphite called graphene rolled up into a cylinder. 
Graphene has a honeycomb lattice formed by covalcntly 
bonded carbon atoms. Its electronic properties are de- 
termined by the p z orbital of the carbon atom. The low 
energy spectrum of graphene consists of two Dirac cones 
located at the K and K' = —K points of the graphene's 
Brillouin zone, where the valence and conduction bands 
touch. To characterize the two Dirac cones, we introduce 
the valley index r — ±1, where t = 1 corresponds to the 
K' point and r = —1 to the K point. The behavior of 
graphene in the presence of an external potential can be 
described by an effective mass approximation, or k ■ p 
theorjES. In this approximation, the envelope wave func- 
tion for the A and B sites of the two-atom unit cell in 
a honeycomb lattice follows an effective Dirac equation. 
In a carbon nanotube, the cylindrical structure imposes 
a quantization condition that leads to either metallic or 
semiconducting nanotubes, depending on the orientation 
of the underlying lattice with respect to the symmetry 
axis of the tube. Here, we will focus on the behavior of 
semiconducting nanotubes. 



A. Single particle spectrum and interactions 
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FIG. 1: (Color online) Schematic representations of a double 
quantum dot in carbon nanotube. Top: Blue regions corre- 
spond to the potential barriers. Bottom: Double dot potential 
at finite detuning A. Dashed lines schematically represent the 
single-particle orbitals. 



In this section we follow previous worlP^EHIHl t derive 
the localized eigenstates of a semiconducting nanotube 
with an additional confinement potential along the tube, 
which is controlled by external gates. We describ e the 
confinement potential of each dot by a square wel l 12 ! 22 ! 
(see Fig. [T]) . The form of the confinement potential will 
not affect our results qualitatively and we note that for 
a single dot, the results for parab olic c onfinement and 
square well are in good agreement! 18 * 20 !. The external 
potential leads to a discrete set of bound states. Taking 
£ as the direction of the nanotube axis and <p as the 
angle perpendicular to the nanotube axis, we can write 
the single particle Hamiltonian as 



H a = -ihvircrx—d^ + a 2 d c ) + V((), 



(1) 



where v is the Fermi velocity, <ii are the Pauli matri- 
ces operating over the sublattice space, and V(() is the 
external potential that describes one or two dots. The 
eigenstates are determined by matching the solutions for 
the dot and barrier regions, which are of the form: 



_ e i(K.R v +kQ 



-K,k 



(2) 



Here k, k denote the wave vectors around and along the 
tube, z T K k — ±(tk — ik)/\/n 2 + k 2 , and the energy is 

given by E K ^ = ±hv^/ k 2 + k 2 . Solving the effective 
Dirac equation for Vm(Oi which is for < C < L 
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and V g otherwise, leads to a quantization condition for 
the longitudinal momentum modes k n of the localized 
states, where n denotes the band index-^l. 

So far, the K and K' solutions are degenerate and 
independent of spin, a =ti, leading to a four-fold sym- 
metry. However, t his symm etry is broken by spin-orbit 
coupling correction d 17 * 23 * 2 ^ and by a constant magnetic 
field B along the nanotube axis £. The spin-orbit cou- 
pling and the presence of a magnetic field modify the 
quantization condition in the (p direction, 

K = k + $ ab /<P R- sA so /(hv) (3) 
k = t/3R. (4) 

Here, s — ±1/2 is the quantum number correspond- 
ing to the spin operator parallel to the nanotube 
(5 c |f) = 1/2 |t) and S c \i) = -1/2 |i)); A so « 
lmeV/i?[nm] is the energy splitting due to spin-orbit cou- 
pling, $ab = BttR 2 is the Aharonov-Bohm flux through 
the nanotube^ and 5>o = hc/\e\. The magnetic field 
also leads to a spin Zeeman term H z = shuj, where 
lj = \e\gB /(2mQc) is the Zeeman frequency in terms of 
the gyromagnetic constant g, the electron mass mo, and 
the speed of light c. The sign convention for the spin 
used here is opposite to the one of Ref. [T5J Note that in 
Eq. Q we only consider the lowest mode in the trans- 
verse direction, since excitations involve energies of about 
hv/R which are much larger than longitudinal-excitation 
energies or Coulomb-interaction effects as long as R -C L. 

An example of the single particle energy spectrum of a 
single dot is shown in Fig. [2] Note that we measure the 
energy with respect to the center of the gap, so that the 
dominant part of the single-particle energy is constant 
and given by Tivk ~ 220 meV/i?[nm]. Generally, disorder 
or the confinement potential itself can lead to intervallcy 
coupling. However, for a noticeable effect, the potential 
must change on the scale of the nearest neighbor lattice 
spacing between carbon atoms do = 1.4 A. The experi- 
ment of Ref. [T7| shows only a tiny valley mixing, and we 
neglect intervalley scattering in this work. 

The longitudinal wave vector depends indirectly on 
spin-valley quantum numbers and on magnetic field and 
these effects are included in the multi band calculations. 
However, since k q 3> $>ab/($ R), Aso/(27ii;), k n , the 
single particle energies can be significantly simplified, 
yielding 

£„, Q « E c n (V g )+E a (5) 
E a = B(j[i orh + 2fi spin s) - A so ts . (6) 

Here we have introduced a single quantum number a = 
(t, s) to describe the spin valley degrees of freedom 
(a = K t,-K I, K' 1). The confinement energy 

E°(V g ) = TivyJ k\ + fc 2 is the single-particle spectrum of 
a dot with potential depth V g ignoring magnetic field, 
spin, and spin-orbit dependences; /i S pm = huj/(2B) — 
h\e\g/(4m c); and fi or i, — {n R\e\ / he) . Within this ap- 
proximation the orbital part of the single-particle states 
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FIG. 2: (Color online) Single particle spectrum as a function 
of B for L = 70 nm, V g = 78 meV, and R = 2.5 nm and 
Aso < 0. Solid curves correspond to E„Kf, dashed curves 
correspond to E n Ki, dash-dotted curves correspond to E nK i^, 
and dotted curves correspond to E nK i^, where n = 1,2,3 is 
the band label. 



separate from the spin- valley part. Figure [2] illustrates 
that, for the parameters studied in this work, spin-valley 
splitting is basically the same for all longitudinal bands, 
and Eq. ([5| provides a good description of the single- 
particle spectrum of a single dot. 

Next, we consider a biased double-dot system, 
schematically presented in Fig. [T] In experiment the 
double quantum dot is formed by applying appropriate 
voltages to external gates and we model the resulting 
confinement potential V2d(C) by a square well poten- 
tial that is -A/2 for -a/2 - L < ( < -a/2, A/2 for 
a/2 < ( < a/2 + L, and V g otherwise. The length of 
the dots is L and a is the width of the interdot barrier. 
As discussed previously we do not expect our results to 
change qualitatively, if a smoother potential is used. In 
the double-dot system at finite detuning, the depths of 
the dots change affecting the single particle energies. In 
the numerical calculation, we determine the eigenspec- 
trum of V2D exactly. The main effect of the detuning 
is an energy shift ±A/2 to the single-particle eigenener- 
gies, where "+" corresponds to right dot, and "— " to the 
left dot. Using Eq. ^ and neglecting interdot tunneling, 
the energies of the localized left and right single-particle 
orbitals are approximately 

« E n(V g ± A/2) + E a « E^Vg) ± A/2 + E a . (7) 

When more than one electron is confined in the single 
or double dot, electron-electron interactions become im- 
portant. The electrons interact through the long-range 
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Coulomb potential 

K(ri,r 2 ) 



(8) 



where kd denotes the dielectric constant. Coulomb inter- 
actions allow for certain off-diagonal matrix elements in 
valley space that are produced by intervalley scattering^!. 
However, these matrix elements are small for quantum 
dots with a size much larger than the interatomic dis- 
tance; they are neglected in this worlP^I 

To obtain an accurate description of interacting few 
electron systems, we extend the single-dot treatment of 
Refs. ll8lT^l to the double dot system. We construct single 
particle orbitals localized in the left and right dots from 
the exact single particle solutions of the double dot and 
then we use them to expand the many-body Hamiltonian 
(see more details in Appendix [A|. The single particle or- 
bitals have a weak dependence on the spin- valley degrees 
of freedom that comes from the dependence of the wave 
vectors k and k on r and s. This leads to a dependence of 
the interaction matrix elements on the spin- valley degrees 
of freedom. However, this dependence is very weak and, 
to a very good approximation, can be neglected. Thus, 
interactions can be considered spin-valley independent 
allowing the separation of the orbital and the spin-valley 
contributions in the two-electron solutions^. 



B. Separating orbital from spin- valley degrees of 
freedom 



Since interactions can be considered diagonal in spin- 
valley degrees of freedom, orbital and spin-valley part 
of the two-electron solutions provide independent contri- 
butions to the energies and the wave functions^. The 
total two-particle wave function must be antisymmetric 
with respect to particle exchange and the symmetry of 
the orbital part must always be opposite to that of the 
spin-valley part. Thus, the two-particle spectrum can 
be grouped according to their orbital symmetry or par- 
ity under particle exchange in multiplets of six states if 
the orbital part is symmetric or ten states if the orbital 
part is antisymmetric. The energy splitting between dif- 
ferent multiplets, called e, depends on the orbital part, 
which is determined by electron-electron interactions and 
longitudinal confinement, and it is generally given by a 
correlated state that is represented as a superposition of 
various two-electron orbital wave functions. The energy 
relations within a multiplet are exclusively determined 
by E a that includes the orbital and spin Zeeman terms 
as well as the spin-orbit coupling [Eq. &\. Therefore, 
the spin-valley part of the wave function always has the 
simple form shown in Table [3j 

The magnetic field dependence of the two multiplets is 
illustrated in Fig. [3j The competition between spin-orbit 
coupling and orbital Zeeman energy leads to a ground 
state crossings in the multiplet with antisymmetric spin- 



-A so 





A SO 


AS states, |tiSi,T2S2) 




\KU K'i)- 
m, JC'f)- 
\K't,K'l)- 




S states, \tiSi,T2S 2 ) + 


\Kl,Kl) 
\KX,K IJ t)+ 
A''t, K'f) 


\K I, K f>+ 
\Kl,K'X) + 
|ift, -ff't> + 
\K'^,K'i) + 


\m,m) 

iff, K'i) + 
\K'i, K'l) 



(a) 



1 




FIG. 3: Table: spin-valley multiplets, here |tisi,t 2 s 2 ) ± = 
(|tisi,t2S2> ± |r 2 s2, tisi))/V2. The states in each multiplet are 
grouped in three columns according to their spin-orbit energy. 
Figures: Schematic magnetic-field dependence of antisymmet- 
ric (a) and symmetric (b) multiplets. 



valley part at a critical magnetic field B c — Aso/2/Wb- 
For B < B c the ground state of the antisymmetric spin- 
valley multiplet [Fig. [3] (a)] is in a superposition of spin 
singlet and triplet, since \K l,K' f} - = i\ KK ') + \V\)~ + 
\KK')~ ||t) + )/v / 2. For B < B c , the lowest state of the 
antisymmetric spin-valley multiplet \K f, K J,) is anti- 
ferromagnetic in spin but ferromagnetic (or polarized) 
in valley space. The lowest state of the symmetric spin- 
valley multiplet \K \., K \) + is ferromagnetic in both spin 
and valley space for all positive magnetic fields. 

The energy difference between the different multiplets e 
and the magnetic field B determines the spin-valley sym- 
metry of the ground state. The thick curves in Figs. [3^ a) 
and (b) correspond to the possible ground states which 
we label according to their spin and valley symmetry. 



C. Model Description 

One of the objectives of this study is to be able to de- 
scribe the evolution of the spectrum as the detuning is 
changed from small to large and the low energy configu- 
rations change from |1X, 1R), i.e. one electron per dot, 
to \2L, OR) or \0L, 272), two electrons in the same dot. 
An accurate description of the strong interaction effects 
requires the inclusion of several single particle bands of 
the double dot system. For example, the behavior of the 
\2L, OR) configurations is expected to be very similar to 
that of the doubly occupied single dot and in that sit- 
uation the strong correlatio ns ne ed to be described by 
many single particle orbitals^^. This situation makes 
the description and the interpretation of the results not 
very intuitive. However, we can significantly simplify the 
description by realizing that no matter how strongly cor- 
related the system is, the parity is a good quantum num- 
ber and the states can be classified according to the state 
parity. Thus, we can model the exact system with a sim- 
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pie effective Hamiltonian in the charge degrees of freedom 
that captures this dependence on parity and describes the 
energetically lowest multiplets of the \1L,1R), \2L,0R) 
and \0L,2R) configurations. 

The charge degrees of freedom of two-electrons in a 
double dot can have three configurations: \2L,0R) ± , 
\IL, 1R) , and |0i, 2/?) ± , where ± characterizes the 
conserved orbital symmetry (i.e. + for antisymmetric 
spin- valley states and - for symmetric spin- valley states) . 
Within this model, interaction effects can be obtained by 
diagonalizing the effective Hamiltonians Hs and Has, 
where S and AS denote symmetric and antisymmetric 
spin-valley state. The complex multiband problem is 
then reduced to simple 3 times three matrices. 



V 



■Vex - A 

-ts 




-t s 
V L r t s 
ts V + V ex 



A 



(9) 



In this limit, energies of the lowest symmetric and anti- 
symmetric spin-valley multiplets are 



E AS 



Esv-A + V- V ex 
E sv -A + V + V ex 



(13) 
(14) 



and 



Thus, the energy splitting e = Eg — E 9 A k, 2V ex in the 
(2,0) configuration is mainly controlled by the exchange 
mechanism. 

This effective model allows a simple and intuitive un- 
derstanding of the underlying physical behavior of the 
double dot system. However, to extract the parameters 
V, Vlr, V ex and t 12 we need to solve exactly the single 
and double dot system. In the next Section, we analyze 
the behavior of the double-dot system at different mag- 
netic field, detuning and interaction regimes by compar- 
ing the exact diagonalization solutions discussed in the 
previous subsection with the model Hamiltonian. From 
this comparison, we extract the parameters of the model. 



Has — 






— tAS 



(10) 



tAS V-Ve, 



These effective Hamiltonians include the onsite and near- 
est neighbor interactions V and Vlr, the tunnelings in 
symmetric ts and antisymmetric tAs configuration and 
the detuning effects A. The dependence of the inter- 
action on the symmetry is introduced by an effective 
exchange term V ex that favors the antisymmetric spin- 
valley configuration. Equations [9] and [10] describe the 
energy related with the orbital part of the wavefunction. 
The total energy also contains the contribution of the 
spin-valley part E S v = J2 a i E i( v g) + E a ]n a , where E a 
was defined in Eq.([6]) and n a denotes the occupation of 
states with spin-valley a. 

To gain qualitative understanding of the interaction 
and tunneling terms in the effective Hamiltonian, we an- 
alyze the limiting behaviors of the low-energy spectrum. 
First, we consider the limit of zero detuning and large 
local Coulomb interactions and we obtain that the low- 
est symmetric and antisymmetric spin-valley multiplets 
have energies: 

E» a Esv + Vls + ..., (11) 

AS sv V-V LR -V ex ' V ' 



E'i 



Esv + Vlr — 



2t 2 



V-V LR + V e , 



(12) 



In the strong interaction regime, our numerical calcula- 
tions indicate that, to a good approximation ts ~ tAs ~ 
ii2 . This approximation allows to obtain a simple expres- 
sion for the energy splitting between different multiplets, 

e = E% 



^t\ 2 V ex /{V - V L r) 2 in the limit of zero 



detuning and for V ex «F - Vlr- 

For a biased double dot system, the single particle en- 
ergies acquire an energy shift of ±A/2 and in the limit 
of large detuning the two electrons occupy the same dot. 



III. RESULTS 

In this section, we analyze the energy spectrum of two 
interacting electrons in a double dot in three different de- 
tuning regimes: (A) small detuning, with one electron in 
each dot; (B) strong detuning, with both electrons in the 
same dot and (C) at the crossover between both regimes. 
In subsection |III D[ we study the transport properties of 
the double dot. 

An important energy scale of the two-particle spec- 
trum is the energy spacing e between the lowest energy 
states with an antisymmetric orbital part and the low- 
est state with a symmetric orbital part at zero magnetic 
field B = 0. In all detuning regimes, we find e > 0, 
in agreement with the Lieb Mattis theorem^. However, 
Coulomb correlations can significantly reduce e. 

In our analysis, we have considered different double 
dot configurations by changing the length and depth of 
the dots as well as the interdot distance and found the 
same scaling of interaction effects as discussed in ReffTSl 
In this section we present results for the parameters 
L = 70 nm, a = 20 nm, V g = 78 meV, R = 2.5 nm. 
A single well with these parameters supports five bound 
states and has an energy splitting between the lowest 
two of them of Tiujq s» 10.6 meV. The dielectric con- 
stant kd is varied between (1.5 <kd< 3.5) which allows 
us to explore the strongly interacting regime where new 
transitions occur. Experimentally, however, it is easier 
to change the length of the dots. The parameter that 
characterize the strength of the interactions is the ratio 
U/(frujo) where U is the characteristic intradot interac- 
tion energy U = e 2 /k^L. This ratio is typically between 
1 < U/(Twjo) < 5 implying a moderate/strong interac- 
tion regime. Another relevant parameter that determines 
the spin-valley nature of the ground state is, as we will 
discuss below, the ratio 2V ex / Ago that reflects the com- 
petition between interactions and spin-orbit effects. 
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Low energy spectrum of symmetric double dot 
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FIG. 5: (Color online) Parameters of the effective Hamilto- 
nian as a function of the dielectric constant. 



FIG. 4: (Color online) Low-energy spectrum of a double dot 
as a function of B. Solid symbols (red online) correspond to 
states belonging to the antisymmetric spin-valley multiplet 
and open symbols (blue online) to the symmetric spin-valley 
multiplet. The solid curves correspond to the effective Hamil- 
tonian description. 



First, we analyze the numerical results obtained with 
the multiband treatment. Figure [4] shows the low-energy 
spectrum as a function of B for kd = 2.5 in a double dot 
with A = (zero detuning) . In the low-energy spectrum, 
we can recognize the symmetric and antisymmetric spin- 
valley multiplets discussed in Table/Fig. [3] The energy 
difference e between the two multiplets is not observable 
in the energy range of Fig. [4] Since interdot tunneling is 
very small compared with the interaction energy, the two 
electrons occupy different dots to avoid strong intradot 
interactions. The interdot interaction is almost indepen- 
dent of the orbitals occupied in each dot. Thus, at small 
detuning there is a negligible occupation of higher bands. 

In Fig. [4j we identify three states that will be rele- 
vant for the discussion of Pauli blockade. One of them 
has mixed valley-spin symmetry and we label it as \M). 
The second state is antiferromagnetic in spin and ferro- 
magnetic in valley degree of freedom and we label it as 
\AfF). The third state, labeled \FF), is ferromagnetic 
in both spin and valley degrees of freedom. States \M), 
\AfF) belong to the antisymmetric spin- valley multiplet 
and \FF ) belongs the symmetric spin- valley multiplet. 
Their configurations are approximately 



\FF) using Eqs. Q, ^TTJ and Q, 



|M) « \lL,lR)+\Kl;K'tY 
\A~fF) « \ILAR)+\K UKty 
\FF) » \1L,1R)-\KI;KI), 



and 



(15) 
(16) 
(17) 



Using exact diagonalization, we extract the wave function 
and conclude that the states li?}^ are given, to a 
very good approximation, by a single Slater determinant 
formed with left and right orbitals in the lowest band. 
We can approximate the energies of \M), \AfF) and 



E 



\M) 



2E1(V g ) - 2t\ 2 /{V - V LR - V e , 



V, 



LR 



(18) 



E, 



\AfF) 



2E c x {V g ) - 2t\ 2 /{V - V LR - V ex ) + V LR 



-2Bfi or b 7 



(19) 



E ]F - F) a 2E c x {V g ) - 2t\ 2 /{V - V LR + V ex ) + V LR 

—Aso — 2B(fi orb + fj, sp in). (20) 

At B — 0, the ground state \M) is only separated by the 
very small superexchange energy e w At\ 2 V ex /{V — Vl R ) 2 
from the lowest spin-valley symmetric triplet that con- 
tains \F F). These four states are separated from the rest 
of the spectrum by a much larger energy scale given by 
Aso- This energy structure resembles the singlet-triplet 
splitting in GaAs double quantum dots. At finite fields, 
\FF) is the ground state. The first excited state changes 
with increasing magnetic field from \M) to \AfF) at 
B c = Aso/(2f^orb), as shown in Fig. [2] This crossing be- 
tween two antisymmetric spin-valley states has no anal- 
ogous in standard GaAs quantum dots. 

From the analysis of the single and double dot spec- 
trum, we can obtain the parameters of the charge ef- 
fective Hamiltonian (Eqs. [9 10). Figure [5] presents the 
parameters V, Vl R , t\i and V ex for a double dot with 
L = 70 nm, a = 20 nm, V g = 78 meV, R = 2.5 nm, 
and A = 0. The V and V ex are obtained from the two- 
electron spectrum in a single dot and Vlr and t\ 2 ob- 
tained from double dot spectrum. The black solid curves 
in Fig. [4] show the prediction from the model using the 
parameters from Fig. [5] 



B. Low energy spectrum for large detuning: Two 
electrons in a single dot 

When the detuning becomes larger than the intradot 
interaction, both electrons occupy the same dot, and the 
charge degree of freedom of the low-energy eigenstates 
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can be described by the \2L, OR) configuration. In this 
regime, the energy spectrum resembles the one obtained 
for two electrons in an isolated dotP. 




FIG. 6: (Color online) Low energy spectrum at large detun- 
ing (A w 35 meV). Solid symbols (red online) correspond to 
states belonging to the antisymmetric spin- valley multiplet, 
and open symbols (blue online) to the symmetric spin-valley 
multiplet. Curves represent the model predictions. 



Figure [6] presents the low energy spectra for large de- 
tuning for kd = 2.5. The solid curves represent the effec- 
tive Hamiltonian predictions. The multiband structure 
of the solutions introduce corrections to the spin-valley 
dependence of the spectrum which can be absorbed in ef- 
fective [i or b and Ago- However, this corrections are small 
and correspond to a few percent changes to the bare /i or h 
and Ago values. 

At zero magnetic field, the antisymmetric spin-valley 
multiplet is favored. This is in agreement with Lieb Mat- 
tis theoreirPH that states that the two particle ground 
state always has a symmetric orbital part. In our two- 
electron system, we can understand this prediction from 
the analysis of the orbital symmetry of the wave function. 
In the noninteracting limit, the orbital ground state is 
constructed with both electrons in the lowest band cor- 
responding to a symmetric orbital wave function, i.e., an 
antisymmetric spin-valley wave function. To form an an- 
tisymmetric orbital wave function at least one electron 
has to occupy an excited state. If the electrons were 
noninteracting, e would be given by the level splitting, 
TujJq, between the first two bands. However, interactions 
substantially reduce the energy difference e between the 
two multiplets as shown in Fig. [6] Consequently, e can 
be changed by tuning the ratio U/Tiu)a, which can be 
achieved by changing the dielectric constant (modifying 
U) or changing the dot length (modifying both U and 
ujq). In the limit of infinite interactions, the electrons are 
strongly localized at the positions that minimize the in- 
teraction energy. Then the orbital symmetry of the wave 
function becomes irrelevant and e vanishes. This effect is 
a signature of a formation of a Wigner molecule^®! 

In the effective Hamiltonian, the formation of a Wigner 
molecule is manifested in the reduction of the parameter 



V ex . This is evident in Fig. [5] that shows that V ex is a 
growing function of the dielectric constant. Even though 
the results of Fig.[5]are for A = 0, we note that V and V ex 
depend weakly on the detuning and can be approximated 
by A = predictions for all detunings studied here A < 
35 meV. In contrast, the tunneling t\i is strongly affected 
by the detuning and is reduced almost by approximately 
a factor of two in comparison with the A = case. 

The reduction of e implies that the exact eigenstates 
become strongly correlated and cannot be written as non- 
interacting wave functions to characterize them. How- 
ever, we can still label the low-lying (2,0) states accord- 
ing to their conserved quantum numbers: 



|M) 
WF) 
\FF) 



\2L,0RY 
\2L,0RY 
\2L,0R)' 



\Kl;K' t>" 

\Ki;Kty 
\Kl;Ki) 



and 



(21) 

(22) 
(23) 



Here the states |2L, CAR)* are correlated orbital states of 
two electrons in the left dot. Our numerical calculations 
shows that several bands are needed to represent these 
states accurately. We note that Secchi and Rontani found 
qualitatively the same correlation effects for a parabolic 
well with weak confinemenlPSl, which suggests the robust- 
ness of these correlation effects. We note that because of 
the conserved symmetries, the eigenstates at small de- 
tuning \M), \AfF), and \FF) will evolve in \M), \AfF), 
and \FF) for large detuning. 

Using the effective Hamiltonian along with the approx- 
imate description of the single particle energies [Eqs(7j 
we can obtain simple expression for the (2,0) 
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configuration energies: 



E m « 2E c 1 (V g ) + V 
E {AfF) « 2E1(V g ) + V 
E\ FF) « 2El(V g ) + V- 



V ex 
V ex 
V ex 



- A so - A, (24) 

- 2Bfj, rb - A, an<l25) 

- A so 

A. (26) 



At B = 0, \M) is the ground state. Above a critical 
magnetic field, there is a ground-state transition to ei- 
ther \AfF) ( Ago < e) or \FF) ( A so > e) due to the 
orbital Zeeman term. Thus, the reduction of e leads to a 
ground state transition to the ferromagnetic state \F F) 
at finite magnetic field. In practice, the formation of a 
ferromagnetic ground state can be experimentally con- 
trolled by changing the length of the dots. 



C. Transition from a double- to a single-dot 
regime. 

We now analyze the transition between the limiting 
behaviors discussed in the previous two sections. 

Figure [7] shows the crossover between \1L,\R) states 
and |2L, OR) states of the double dot system at B = 0. 
In absence of interdot tunneling, the two-electron spec- 
trum shows sharp crossings between the \1L, 1R) states 




23 24 25 26 27 28 

A [meV] 



FIG. 7: (Color online) Zoom of first crossing of (1,1) and (2,0) 
states for B = 0. Solid symbols (red online) correspond to 
states belonging to the antisymmetric spin- valley multiplet, 
and open symbols (blue online) to the symmetric spin-valley 
multiplet. Curves represent the model predictions. 



and |2L,0i?) states with increasing detuning. Because 
of the interdot tunneling, crossings between states with 
the same symmetries turn into avoided crossings. The 
avoided crossings occur at the same critical detuning 
for states of the same multiplet. The avoided crossings 
within the multiplet with an antisymmetric orbital part 
occur at relatively larger detuning since the tunneling 
electron is forced to occupy an excited band. 

The crossover regime, presented in Fig. [7j is strongly 
affected by interactions. The difference between the criti- 
cal detunings belonging to the avoided crossings of states 
with symmetric and antisymmetric orbital part is a di- 
rect measure for the energy splitting e. Furthermore, 
correlations in the \2L, OR) states decrease the tunneling 
coupling to the corresponding \1L, 1R) state, leading to 
a sharper avoided crossing. 




23 24 25 26 27 
A [meV] 



FIG. 8: (Color online) e as a function of detuning, symbols 
correspond to the numerical results and the curve is the ef- 
fective Hamiltonian prediction. 



In the effective Hamiltonian description, the lowest 
\1L, 1R) and |2L,0i?) configurations are close in energy 
when A ~ V — Vlr- In this regime, the \0L, 2R) con- 
figurations are energetically suppressed and do not affect 
the low energy spectrum. Within this approximation, we 
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can obtain a simple expression for e, 

e = V ex + i (\J 4<J 2 + (A — V ex — V + Vlr) 2 

-^U\ 2 + (A + V ex - V + V LR )^j (27) 

This expression compares well with the numerical pre- 
diction as shown in Fig. [8| 




0.0 0.1 0.2 0.3 0.4 0.0 0.1 0.2 0.3 0.4 



B [T] B [T] 

FIG. 9: (Color online) Phase diagram for a double quantum 
dot with five bands and kd — 1.5 (a) and kd = 3.5 (b). Black 
solid lines separate regions with different spin- valley symme- 
try. Regions with the same greyscale ( color online) charac- 
terize the ground state. 



From the analysis of the low-energy spectrum, e. g., 
Fig. [7J we can extract the phase diagram as a function 
of the detuning A and the magnetic field B. The phase 
diagrams for strong (kd — 1.5) and weak (kd — 3.5) inter- 
actions, presented in Fig. |9j exhibit clear differences. At 
exactly zero magnetic field, the ground state always has 
a symmetric orbital wave function (\M) and |M)). This 
state is preferred with respect to other states within the 
same multiplet by the spin-orbit coupling Ago and, with 
respect to states of the other multiplet, by the energy gap 
e. The ground state changes at a critical magnetic field 
to a valley-polarized state. The valley polarized ground 
state has symmetric orbital part for Ago < e and an an- 
tisymmetric orbital part for Ago > e. For small detuning 
A, the energy splitting e is caused by the small superex- 
change and e < Ago- However, for large detuning, both 
electrons are on the same dot, and the value of e can 
be larger or smaller than Ago. Thus, if 2V ex > A$o, 
there is a transition to a valley-polarized state with a 
symmetric orbital part, \AfF), as detuning is increased 
[as observed in Fig. |j] (b)]. 



D. Sequential transport 

In this section, we discuss how the two-electron cigen- 
spectrum affects sequential transport through a double 
dot. In particular, we analyze the existence of Pauli 
blockade in serially coupled dots that would lead to a 
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current rectification in DC transport!^. Figure 10 shows 
a schematic description of two situations when a Pauli 
blockade occurs. In both cases, the double-dot system is 
detuned such that there is always at least one electron 
in the left dot. States relevant for transport through the 
double dot are the \2L,0R) states \M) and \FF) as well 
as the \1L, 1R) states \M) and \FF). In Fig.fioj the ver- 



tical axes denotes energy and the horizontal axes denotes 
the spatial coordinate along the nanotube. The center of 
each figure shows the double dot created by tunnel bar- 
riers to the contacts and between the dots. The gray 
rectangles to the left and right of the double dot are the 
Fermi seas in the contacts. 



(a)B=0 Blocked 



~4r ~ 



\M) 



\F~F) 



\M) 



(b) B=B* Blocked 

|M)n 



|FF> 



m 



|FF) 



FIG. 10: (Color online) Schematic representation of a Pauli 
blockade for: (a) zero magnetic field and negative bias, and 
(b) finite magnetic field and positive bias. The white areas 
in the center of each figure correspond to the two dots. The 
black circle in the left dot stands for the position and energy 
of one of the electrons. The states \M) and \M) are repre- 
sented by solid lines, and the states \FF) and \FF) by dashed 
lines. These states are organized vertically according to their 
energies, and also represent the possible positions of the sec- 
ond electron. The wide gray rectangles next to the left and 
right dots correspond to the energy bias. 



For a positive bias voltage (opposite to Fig. 10 (a)) , the 
electrochemical potential in the left contact is larger than 
in the right one, and current flows via a sequential tunnel- 
ing process \1L,0R) ->■ \2L,0R) -> \1L,1R) -> \1L,0R). 
Interdot tunneling is assumed to conserve the spin-valley 
degree of freedom and allows for transitions between 
states \M) and \M) or \FF) and \FF). Because of the se- 
quential transport setup, the left dot couples to the left 
reservoir, allowing for transitions between the |1Z, OR) 
and \2L,0R) states. Analogously, the right contact al- 
lows for transitions between the \1L, OR) and \1L, 1R) 
states. Figure 10 (a) considers Pauli blockade at B = 0. 

In this scenario, the current is blocked when the state 
\FF) is occupied by an electron tunneling in from the 
right reservoir. Once in state \FF), the electron cannot 
tunnel back because of the filled Fermi sea in the right 
contact. Also, if e is large, the electron cannot tunnel 
to the left dot since a transition to the \FF) state is 
energetically suppressed [see Fig. 10 (a)]. However, e can 
be strongly reduced, allowing for a finite exit rate from 
state \FF). This effect might explain the absence of a 
Pauli blockade in Ref. [TTJ This Pauli blockade implies 



a rectification of the current since, by inverting the bias 
voltage, a finite current can fkmW. 

In order to make these statements more quantitatively 
we calculate the stationary current with a rate equa- 
tion approacli22H3U. We describe the regime of possible 
blockade depicted in Fig 10 (a) while assuming spin-orbit 
coupling to be much bigger than temperature, transport 
voltage and external coupling the the reservoirs. 

The rate equations and the resulting stationary current 
are given in Appendix [B] The dependence of the station- 
ary current Iu in the blockade setup and the probability 
to be in any of the three degenerate \FF) states are given 
by: 



FF 



B + C{2t 2 +e 2 ) 
C{t 2 +e 2 ) + D 
B + C{2t 2 +e 2 ) 



(28) 
(29) 



Here e denotes the energy splitting between the states 
\M) and \FF) and tg is the interdot tunneling rate be- 
tween states \FF) and \FF). The constants A,B,C, D 
depend on the coupling strength to the contacts and spec- 
tral weights for the tunneling probabilities and their form 
is given in Appendix |B| 

If e is the dominant energy scale then the current is 
suppressed as Iu oc 1/e 2 and the double dot gets stuck 
in states \FF) with P F ~ F — > 1. This is the regime of 
Pauli-blockade. However, we find that due to interaction 
effects interdot tunneling t$ can even exceed the energy 
splitting e, thus removing the blockade mechanism. 




e [meV] 



FIG. 11: (Color online) Ratio between leakage current Iu ( 28 I 
corresponding to setup in Fig. 10 (a) and the finite current 



for reversed transport voltage I op |B3j), i.e., the current in 
the open (unblocked) direction. Used parameters Tl — T_r = 
O.OlmeV, t A s =t s = O.lmeV and S M = 0.7 S F = 1. 



Due to Pauli-blockade the double dot acts as a cur- 
rent rectifier, since by reversing the transport voltage in 
Fig. [To] a) a finite current can flow. 



Fig. 
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illustrates 

how effective the double dot acts as current rectifier as 
function of e. 

At finite magnetic field, the transport behavior can 
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also be very interesting. For example, if interactions are 
strong enough to suppress e to become on the order of the 
spin-orbit coupling (A so > e) but still much larger than 
the interdot tunneling (e 3> t), then a current blockade 
can occur in the opposite bias direction of the zero field 
case. This situation, depicted in Fig. 10 (b), is achieved 
by changing both the magnetic field and the detuning. A 
magnetic field B* e/2(/i or fc + (J, s pi n ) is applied so that 
states \FF) and \M) are degenerate, while in the |1L, 1R) 
charge configuration, the state \FF) is ground state and 
separated from state \M) by the energy e (because of 
the applied field B*). The physical situation of Fig. 10 
(b) resembles that of Fig. 10 (a), and also results in a 



current blockade. This ability to control the direction 
of the current rectification by varying the magnetic field 
and detuning can have applications to carbon nanotube- 
based spintronic proposals. 



E. ID-Quantum dot arrays 

The rich physics associated with the strong correla- 
tions of double occupied dots is expected to have sig- 
nificant impact in the many-body behavior of a one- 
dimensional chain of dots. The double dot analysis car- 
ried out here can be used as a starting point to study 
the behavior of a linear array of coupled dots. Accord- 
ing to our analysis, there exists an interaction regime in 
which the behavior of two electrons in the same dot is 
strongly correlated and controlled by e while the behav- 
ior of electrons in different dots is weakly correlated and 
can be accurately described considering only the lowest 
orbital of each dot. In order to generate an effective 
Hamiltonian that captures the essential effects of strong 
onsite correlations, we apply the ideas of Hubbard oper- 
ators to describe double occupied dots or doublons (see 
i.e. Ref. This description assumes that the strong 
onsite Coulomb repulsion suppresses the probability to 
occupy a dot with more than two electrons. This is a 
good approximation for low enough fillings. 

Formally we start from the complete Hamiltonian that 
describes a chain of dots. This Hamiltonian represents 
an extension of the double dot Hamiltonian introduced 
in the Appendix. 

We introduce the operators a\ n a which create an elec- 
tron in the r dot, in the n orbital state and with a given 
spin-valley configuration a. To simplify the notation we 
introduce a], a = a\. 1 a for electrons in the lowest band 
which describe well the singly occupied dots. The state 
that describes two electrons in the r dot can be expanded 
in single particle basis as 



\d>r),r,a,a') — f a ,a' flr],n,m { a ]-,n,a a \,rn,a' 

+ (-l) r '«La'< m ,J|0) (30) 

Here, rj labels the spin-valley symmetry of the two- 



electron state r] = 1 for antisymmetric spin-valley states 
and r] — 2 for symmetric spin-valley states and f a .a' = 

defines the doublon opera- 

r,a, a >\°) and 4,,r, a , a > = 

3r],n,m do not depend on 
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yl + <5a,a'- Equation 
tor such that {d^^^.a' , 

(- 1 )' 7 4,T-,a',a- Note that the I 

a and a'. The coefficients l3 Vt n tm can be obtained by 
diagonalizing the local part of the Hamiltonian which 
amounts to solve two electrons in a single dot. We now 
assume that the occupation of higher excited two-particle 
states, which are separated by an energy gap of the or- 
der of the single particle level spacing, can be neglected. 
This should be a good approximation for the low energy 
spectrum. 

The effective Hamiltonian is obtained by projecting the 
complete Hamiltonian onto the subspace of empty, singly 
and doubly occupied dots and can be written as 

H ef f =P(H e + H d + H ed + V C )P. (31) 

where P = Yii Pi represent a projector to the physi- 
cally allowed subspace where Pi = |0)i(0| +J2 a \ a )i{ a \ + 
^2 \d rj )i(d rj \ (see Ref. l33|) . H e describes the behavior of 
single occupied dots, H d describes the behavior of dou- 
ble occupied, H ed contains a coupling between singly and 
doubly occupied dots and V c represents the Coulomb in- 
teraction between different sites. The explicit form of 
these contribution is: 

H e = ^ F r a o) ra (l ra -\- ^ tr,r' &r a ®'r' a (32) 
r,a (r,r'),a 

(r,r') . 77, a<a' 

H e d — ~y K ^yf aiCl3 f a2a3 d\j raia3 a} r i a2 a rai d^ r i a2a3 

(r,r) ,rj,rj' , 
ai,a2,Q3 

~l~ ^ fact' 9rj^r.r'\d^ r aa iQ ir a Q ir f a t -\- H.c] 

{r,r')r],a : oi' 



(34) 



The diagonal part of H e and H d represent the elec- 
tron and doublons energies composed of the spin-valley 
energy and onsite energies {E r ^ a = Ef(V r )+E a and 

<«V =2EZ(V r )+E a + E a ,+V+(-l)"V ex ) , andtand 
v- d > the tunnelings of the electrons and doublons. The 
H ed contains a term that describes the hopping of an elec- 
tron from the doubly occupied site to a singly occupied 
site which in this language corresponds to an exchange 
of a doublon and an electron. The second term in H ed 
represents the hopping of an electron from a singly occu- 
pied site to a singly occupied site and viceversa, which is 
represented as destroying two atoms and creating a dou- 
blon. Finally, we have the offsite Coulomb interaction V c 
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in terms of the dot density n r denned as 



E 



a ru a ru 



'> E 

r),{a,a' } 



(f d 

r},raa' r),raa 



(35) 



Here, the summation {a, a'} is restricted to a < a' for 
rj = 1 and a < a' for rj = 2. For the offsitc interaction, 
we assume exclusively capacitive coupling and neglect a 
dependence on the symmetry of the two particle states. 

Explicit expressions for the parameters g v ,r,r'j and t^, 
can be obtained by comparing matrix elements of the ex- 
act and effective Hami hton ians If the doublon solution 
is expressed as in Eq. 
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the g n ,r,r'i an d t^, can be 
expanded in terms of the"" doublon expansion coefficient 
Pri,n,m and the many-body Hamiltonian matrix elements 
(see e.g. Ref. [33]). Alternatively, these parameters can 
be obtained from comparison between exact and effec- 
tive Hamiltonian solutions for two and three-electrons 
systems. These few-electron calculations might be chal- 
lenging but allow the determination of the parameters 
needed for many-body calculations. The extraction of 
these coupling parameters is beyond the scope of the cur- 
rent article. 

In the most general case the parameters that describe 
the effective Hamiltonian depend on lattice sites posi- 
tions and can be controlled by changing the detuning in 
each lattice site A. r . For example, an enhancement of the 
superexchange interactions can be achieved by detuning 
some of the dots and, therefore, reducing the energy cost 
of double occupancies. 

The effective Hamiltonian can be applied to describe 
an array of coupled dots in many different regimes and 
its phase diagram might exhibit novel phases. In par- 
ticular, the existence of spin-valley "triplet" close in en- 
ergy to the spin- valley "singlet" can lead to phenomena 
richer in comparison to the standard single-band Hub- 
bard modeP^. For example, for one particle per site (fill- 
ing n = 1), it is known that the ground state of the 
usual Hubbard model has infinite susceptibility to spin 
dimerizatiorPS i.e. formation of singlet /triplet bonds be- 
tween nearest-neighbor sites. However, the ground state 
is not dimerized since the formation of singlet bonds, say 
at sites (2 i, 2 i+1), is penalized by the large energy cost 
of the remainder triplet components between sites (2 i+1, 
2 i+2). In carbon nanotubes operating in the strongly in- 
teracting regime , e « 1, even spin- valley "triplet" states 
can lower their energy by virtual hoping. The later sit- 
uation reduces the energy cost of triplet formation and 
the infinite susceptibilities to spin dimcrization might in 
this case translate in actual dimerization of the ground 
state. 

Away from n = 1, the presence of a spin- valley 
"triplet" can also significantly affect the magnetic struc- 
ture of the system. In particular, for systems which 
already exhibit ferromagnetism in the standard single- 
band Hubbard model, the inclusion of the spin-valley 
"triplet" can strengthen and extend the ferromagnetic 
phase. For example, generic 2 and 3 dimensional square 



lattice geometries, and others with similar connectivity 
conditions, exhibit Nagaoka ferromagnetisnpSI when the 
gain of kinetic energy of a single hole exceeds the decrease 
of superexchange energy. The latter condition is fulfilled 
at very large on-site interactions. In quantum dots in 
carbon nanotubes with small e, Nagaoka-type ferromag- 
netism might become stable at reduced value of the inter- 
action since the virtual hopping of the spin- valley triplets 
reduces the super-exchange penalty of having a polarized 
state. 

For arrays with filling 1 < n < 2 and zero detuning, the 
low energy physics consists of doubly and singly occupied 
dots while the occurrence of empty sites is strongly sup- 
pressed since it implies an increase of double occupancies. 
This implies that the number of electrons and doublons 
will be independently conserved and that the second term 
of H e( i can be neglected. Also, the tunneling of the elec- 
trons and the doublons can be neglected. Thus, the only 
relevant terms in the effective Hamiltonian are the on- 
site energies, the long-range Coulomb and the electron- 
doublon exchange (first term in H ec i). This regime can 
lead to interesting phenomena when the electron-doublon 
exchange becomes comparable the multiplet splitting e. 

Finally, it should be pointed out that the long-range 
Coulomb interaction can have a preponderant influence 
in the charge distribution in small arrays of quantum 
dots. 



IV. SUMMARY AND CONCLUSIONS 

We have presented a detailed study of the few-electron 
eigenspectrum of a double quantum dot in a semicon- 
ducting carbon nanotube. We showed how the spin- 
valley physics leads to the formation of multiplets. The 
internal energy structure of the multiplet is practically 
unaffected by a change in either the confinement poten- 
tial or the interaction strength, but the energy gap be- 
tween different multiplets is strongly modified by both. 
We showed that for sufficiently strong interactions, the 
spin-orbit coupling can exceed the energy splitting be- 
tween states with symmetric and antisymmetric orbital 
parts for any detuning between the dots. This situa- 
tion modifies the two particle phase diagram. Above a 
critical interaction strength, the ground state at small, fi- 
nite magnetic fields is always ferromagnetic independent 
of the detuning. Furthermore, in this strongly interact- 
ing regime, the blockade of linear transport gradually 
disappears since the reduction of the multiplets' energy 
splitting e allows a finite tunneling probability from the 
\1L, 1R) state to the \2L,0R) state even for a symmet- 
ric spin- valley part. Pauli-blockade physics will occur for 
weak enough Coulomb correlations, which can be sup- 
pressed by either working with short dots, or by covering 
the nanotube by strong dielectrics^. We note that a well 
developed Pauli-blockade is the precondition for realiza- 
tion of spin-qubits in double dots allowing for coherent 
singlet-triplet manipulation^^!. 
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Our understanding of the double dot physics can be 
used as an starting point to analyze the behavior of an 
array of coupled dots. The effective Hamiltonian, as pre- 
sented in Eq. [3TJ can be applied to describe an array of 
coupled dots in certain regimes. In the future, we would 
like to explore the many-body physics associated with 
the strong onsite correlations and its consequences in the 
magnetic ordering of the system. 

In the analysis presented here, we neglected terms that 
flip either the spin or valley degree of freedom, or both. 
However, such mechanisms could change the various level 
crossings to avoided crossings and thus open new ways to 
control the spin-valley degree of freedom as well as the 
transport properties of this systenPIH13 
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Appendix A: Coulomb matrix elements 



To evaluate the interaction matrix elements, we use 
the following relation: 



1 _ 1 

r 7r 



(A4) 

If the rest of the integrand does not depend on <p\ and 
if 2, the integration over those coordinates can be done 
analytically, 



1 

4^ 



27T r 2n 



K (i?|2qsin[(<pi - ip 2 )/2]\)dip 1 dip 2 = 

l (R\q\)K (R\q\). (A5) 



Also, the integration over d and £2 can be done ana- 
lytically since $(C) has a simple functional dependence 
that can be written in terms of exponentials with real 
or imaginary arguments. Once those integrals are ob- 
tained, we are left with the momentum integration over 
q in Eq. (A4). Defining 



F np (q) 



- / JiCi 



1 {*;(Ci)<MCi)Ri, (A6) 



we can write the Coulomb interaction matrix as a one- 
dimensional numerical integration over the momentum 

q, 



We create an orthonormal set of localized single parti- 
cle orbitals <J>i solving exactly the double dot system and 
we use this set to construct the many-body Hamiltonian 
for the interacting double-dot system: 



i¥=3 



a k ai. 



ijkl 



(Al) 



Here Ei = (<&j|i?o|5?i) are the single-particle energies, 
Tij = ($>i\Ho\$j) are the tunneling matrix elements, 
and Uijki — ($i$j\V c \$k&i) correspond to the interac- 
tion matrix elements. At zero detuning, the tunneling 
matrix elements, Ty, are only nonzero between states 
with the same a, r, n and located at different wells. For 
nonzero detuning, tunneling between different bands be- 
come possible. The electrons interact through the long- 
range Coulomb potential, Eq. [HJ and the interaction ma- 
trix elements take the form 



U, 



nrapr 



All 2 k. 



- I <(Ci)-* P (Ci) 

d J 



C(C 2 )-$r(C2)-#ldClWC2, (A2) 

r 



where 
1 1 



r |r x -r 2 | ^2i? 2 [l - cos(^ - V2 )] + (Ci - C2) 2 ' 

(A3) 



U, 



nmpr 



k d TT 



F np (q)F mr (-q) l (R\q\)K (R\q\)dq. 



In the calculation of the Coulomb matrix elements, we 
neglect the dependence of the longitudinal wavevec- 
tor k n on spin, valley, and magnetic field. As dis- 
cussed in Sec. II A| this approximation is good for k ^> 
&ab/(&oR), Aso/(27m>), which is certainly true for the 
parameters used in this paper. 

There is another effective interaction that takes into 
account intervalley scattering: 



U., 



KK' 



1 



<(Ci)-<MCi) 



nmpr ^ 2 

*L(C 2 ) • $ r (C 2 )^'^idCi^ 2 dC 2 (A8) 

where V KK is a short range interaction that can be taken 
proportional to <5(^i — (,2)- We have estimate the U££ 
terms and we have concluded that they only introduces 
minor corrections. For that reason, we neglect this con- 
tribution in the present work. 



Appendix B: Rate equations 

Here we derive the rate equations and the stationary 
current corresponding the the setup depicted in Figfi0| 
a). At zero magnetic field and large detuning the sin- 
gle particle ground state consists of the Kramer doublet 
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\L1K' t) and \L1K |) localized in the left dot. The two- 
particle states are characterized by their charge degree of 
freedom and by their spin valley symmetry The \2L, OR) 
and \1L,1R) states with antisymmetric spin- valley part 
are non-degenarate and are called \M) and \M) respec- 
tively. They are defined in Eq.(21| and ( |15[ ). According 
to Figj3 the \2L, OR) and \1L, liZ) states with symmetric 
spin-valley part are each three fold degenerate and in- 
clude the states \FF) and \FF) respectively. Following 
Eqs( 17 ) and ( 23 ) the structure of these states is given by: 



Pf = P, 



■Pn 



-\F + )^-r\F )^P\F-,, r _ r + 

The corresponding rate equations are given by 



. pF _ pF+ 



Pf°+P^:, etc. 



F-) = 


\FF) = 


2L,0R)- \K l;K I) 


\Fo) = 


\2L,0R) 


-\Ki-K 1 t) + 


F + ) = 


\2L,0R) 


~\K' f,K' t) 


F-) = 


\FF) = 


1L,IR)- \K i;K i) 


\Fo) = 


\1L,1R) 


~\KUK' t> + 


F + ) = 


\IL,1R) 


~\K' f,K' t> 



Due to the serial setup each contact couples only to 
its adjacent dot. The tunneling rates to the collector 
(left contact in FigjlOja)) depend on the correlated two 
particle states, which is accounted for by introducing the 
spectral weights S F and Sm that determine how easy 
an electron tunneling between the left dot and the left 
reservoir can cause a transition between a two particle 
and one particle state on the dot. 

Sp = J2\( L1K i\ a L m Kl\F-)\ 2 
rn 

S M = J2\( L1K i\aLmK't\M)\ 2 



Due to the simple form of the \1L, 1R) wavefunctions 
the spectral weights corresponding to the tunneling be- 
tween right contact and right dot has trivial spectral 
weights S a = I {L1K i\ a mK > t \M) | 2 = \ and S F = 
| (L1K i\ a mK i \FJ) | 2 = 1. 

We denote the reduced density matrix for the dou- 
ble dot by P and denote off-diagonal matrix elements 
with :— (s'\P\s) = (F s s )*, where s and s' are eigen- 
states of the double dot and diagonal elements with 
P s := (s\P\s). We neglect pr ocesse s that can mix states 
with different spin- valley part! 39 * 41 ! since hyperfine inter- 
action is small in carbon nanotube and spin-orbit cou- 
pling is mainly a spin-valley coupling with a very weak 
effect on the longitudinal wave function. Then degen- 
erate states are equally populated and the relevant ma- 
trix elements of the reduced density matrix for the dou- 
ble dot system are given by Pi = P\lxk>\) + P\lik\)\ 



d t Pi 
d t P F 



d t P F 



d t P^ 
d t P M 
dtP M 

J pM 



= -2T e P 1 +T c S f P f + 2T c S m Pm (Bl) 

= -T c S F P F ~it s (P^ -Pf) 

= \T E P 1 + l t s {P?-Pl) 

- {-\v c S F +ie)Pi -it s {P F -P F ) 

= -2T c S M P M -it AS {P^-P§) 

= -r c s M p£ -it AS {P M -p M ) 



with P 1 +P F + P F + P M + P M = 1. Here T E denotes the 
coupling strength to the emitter contact (right contact 
in Fig 10 a) and Tq denotes the coupling strength to the 
collector/left contact. 

These rate equations are solved analytically. The re- 
sult for the stationary current and the occupation of the 
states are given by: 



hi = 



P - — 

FF ~ 

N = 



2 eT ^ Pl = 



l&S F S M t 2 s t 2 AS T c /N (B2) 
T 2 n S F +4e 2 )/N 



h n 

12S Mt 2 AS^ E{2t 2 g + £ 2 ) + 3S F SMtAS^C^E 



2 r 2 
AS L C 
2 i - i q2 t^2 



-2S F t S (AS as^ c + ^as^e + S M T C T E ) 



Here tAs denotes the interdot tunneling between states 
\M) and \M) and ig the interdot tunneling between the 
the states \FF) and \FF). Spectral weights S F and Sm 
account for the overlap between an electron hopping onto 
the dot and the two particle eigenstates, and Tc and T E 
are the coupling strength to emitter and collector. 



The stationary current (B2) vanishes as 1/e 2 for large 
e. For reversed transport voltage, still a finite current can 
pass through the double dot which implies that, for large 
e, the double dot acts as a current rectifier. We calculate 
the current for the unblocked voltage direction. Inverting 
the bias voltage depicted in Fig 10 a), the current flows 



by electrons hopping from the left reservoir through the 
double dot to the right reservoir. Applying again the rate 
equation approach, we obtain: 



hp — 



-P 



2eTit 



2 j.2 
R L AS 



h2t 2 AS (T R + 2S M r L ) + s M r L r 



£3) 



We note that Tl corresponds to Tc in Eq.(B2) and Tr 
corresponds to T E . 
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